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Field of the Invention 

The present invention relates to CDMA (Code Division 
Multiple Access) spread spectrum receivers, and more 
particularly to fast acquisition GPS (Global Positioning System) 
receivers (or, more generally, ranging receivers) . 

Background of the Invention 

Spread spectrum communication in its basic form is a method 
of taking a data signal that is used to modulate a sinusoidal 
carrier and then spreading its bandwidth to a much larger value, 
e.g. in a global positioning system (GPS) application, the 
spreading is achieved by multiplying a single- frequency carrier 
by a high-rate binary (-1,1) pseudo-random noise (PRN) code 
sequence known to GPS users. Thus, the signal that is 
transmitted includes a data component, a PRN component, and a 
(sinusoidal) carrier component. 

At the receiver of such a signal, a synchronized replica of 
the transmitted PRN code is required to de- spread the data 
sequence. Initial synchronization, called acquisition, is 
followed by fine synchronization, which is called tracking. 

The present invention relates primarily to acquisition, not 
tracking. Acquisition is the process by which the replica PRN 
code is synchronized (to within a small timing offset) with the 
code conveyed by the received signal either for the first time 
or after losing a previously acquired signal, and also by which 
the carrier frequency of the received signal is determined. 



Thus, to acquire a signal, an acquisition system must accurately 
determine any frequency- shifting of the received signal from the 
transmitted frequency in order to accurately wipe off (remove) 
the carrier signal. Frequency-shifting can be caused by 
5 relative motion of the transmitter and receiver (Doppler- 

shifting) as well as by clock inaccuracies (so that a 
transmitter and receiver sometimes do not agree on what is in 
fact the same frequency) . The carrier frequency- shifting 
results in a modulation of a code component after carrier wipe- 
10 off in the receiver. Thus, in acquiring a signal, it is also 

necessary that the replica code sequence be not only time- 

Q aligned with the received code sequence, but also modulated to 
compensate for the frequency- shifting so as to fully eliminate 

fl the PRN sequence and leave behind only the data conveyed by the 
15t;' received signal. The acquisition process is therefore a two- 

^|;| dimensional search, a search both in code phase and in 

: . frequency . 

j:55;:i; 

For some GPS signals (the civil GPS or C/A (course 
ni acquisition) code signals) , the search interval in the frequency 

20^^ domain can be as large as +/-6 kHz. In addition, the phase of 
the received code relative to the replica can be any possible 
value of code phase, due to uncertainties in position of the 
satellite and time of transmission of the received signal. A 
PRN code period is typically 1023 chips, the term chips being 

25 used to designate bits of code conveyed by the transmitted 
signal, as opposed to bits of data. Thus, the acquisition 
module of a receiver must search a 12 kHz -wide interval with 

1023xks different code phases, where kg denotes the number of 
samples per chip. 

30 Ordinary GPS receivers, i.e. those designed only for 

operation with unobstructed satellites, search for the frequency 
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shift with a granularity of around 1 kHz. Thus, such receivers 
must search 12xl023xks different code/f requency combinations. 

A GPS receiver designed for indoor operation must have an 
operating mode with equivalent noise bandwidth on the order of 
5 10 Hz in the acquisition stage. Even with an equivalent noise 

bandwidth as small as 10 Hz though, for reliable tracking some 
post -detection filtering must still be performed as well as some 
further refining of the ^value determined for the carrier 
frequency in the acquisition stage. The granularity of 10 Hz 
10 requires that the receiver search 1200xksXl023 different 

code/f requency combinations and makes the sequential search so 



O time-consuming as to be unrealistic, motivating the use of 
ff'fe parallel and fast search methods. 

%s? r 

J{ Because of the granularity used for the carrier frequency 

I5CI in the acquisition stage, the tracking stage often includes a 
1^. preliminary fine frequency process for refining the carrier 
frequency determination, before initiating tracking. An 

alternative is to provide a more precise signal acquisition, but 
to do so is problematic, because signal acquisition is based on 



20' performing a correlation of the received signal with a replica, 

a correlation that includes an estimate of the offset of the 
carrier frequency from a nominal carrier frequency (due to the 
Doppler shifting and clock inaccuracies mentioned above) , and 
since the received signal is modulated not only by the PRN code 

25 (for example at 1.023 MHz for a coarse acquisition signal) but 

also by the navigation message at a bit rate of 50Hz, the 
correlation becomes corrupted if a signal fragment is used that 
is longer than 20 ms (the duration of a navigation data bit) . 
Thus there is an apparent limit of 50 Hz for the precision with 

30 which the carrier frequency can be acquired. It is possible to 
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overcome the apparent limit by performing multiple correlations, 
but then signal acquisition takes longer. 

What is needed is a method of fast acquisition by a GPS 
receiver (or, more generally, any ranging receiver, i.e. not 
necessarily a ranging receiver used with the Global Positioning 
System but including for example ranging receivers used with 
GLONASS, the Russian version of a global positioning system) 
that includes a more precise estimate of the carrier frequency 
than is typically provided (i.e. to within a few Hz, as opposed 
to 50Hz) , and so abbreviate or eliminate the need for any 
preliminary fine frequency process in the tracking loop of the 
receiver. 

Summary Of The Invention 

Accordingly, in a first aspect of the invention, a method 
is provided for determining information about the carrier 
frequency of a signal transmitted by a possibly moving 
transmitter, the signal having a code component and a carrier 
component, the method including: a step of responding to 
successive approximately carrier-demodulated received signal 
fragments, and providing a set of correlation results indicating 
information about the correlation of the approximately carrier- 
demodulated received signal fragments with a replica of the code 
component and any remaining carrier component, wherein the set 
is formed using different possible offsets from a nominal 
carrier frequency used to approximately carrier-demodulate the 
received signal fragment, and further wherein each element of 
the set is provided as a phasor having a magnitude and a phase; 
and a step of responding to the set of phasors, selecting the 
phasor having a magnitude distinguishing it from all the other 



elements of the set, and determining the phase of the selected 
phasor . 

Further in accord with the first aspect of the invention, 
the set of correlation results may be a matrix of correlation 
5 results, and the matrix of correlation results may be spanned by 
an index indicating an offset from a nominal carrier frequency 
and also by an index indicating code phase, and the selected 
phasor may be the phasor having the maximum magnitude of all the 
elements of the set. Further, the step of providing the matrix 
10 of correlation results may include a step of performing a 
H coherent integration of each of a series of signal fragments, 
□ and a step of performing a non-coherent integration in which the 
^- phasor results of the coherent integrations are combined without 
j:; regard to phase. Further still, the step of performing the non- 

157: coherent integration may involve multiplying each element of a 

'if 

SI matrix of correlation results provided using a coherent 

integration of a first signal fragment, by the complex conjugate 
of a corresponding element for an immediately preceding signal 
pj fragment. Also further, in providing the matrix of correlation 

2 OH results as phasor values and in determining the phase of the 

phasor having the maximum magnitude of all the elements of the 
matrix, in some applications only at most two phasor values may 
be held in a memory device at any instant of time, and of the 
two phasor values, only the phasor value having the larger 

25 magnitude may be saved in the memory device before calculating a 

next phasor value. 

In a second aspect of the invention, an apparatus is 
provided for determining information about the carrier frequency 
of a signal transmitted by a possibly moving transmitter, the 
3 0 signal having a code component and a carrier component, the 

apparatus including: means, responsive to approximately carrier- 
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demodulated received signal fragments , for providing a set of 
correlation results indicating information about the correlation 
of the approximately carrier-demodulated received signal 
fragments with a replica of the code component and any remaining 
carrier component, wherein the set is formed using different 
possible offsets from a nominal carrier frequency used to 
approximately carrier-demodulate the received signal fragment, 
and further wherein each element of the set is provided as a 
phasor having a phase and a magnitude; and means, responsive to 
the set of phasors, for selecting the phasor having a magnitude 
distinguishing it from all the other elements of the set, and 
determining the phase of the selected phasor, and for providing 
information about the carrier frequency based on the phase of 
the selected phasor. 

Further in accord with the second aspect of the invention, 
the set of correlation results (phasors) may be a matrix of 
correlation results, and further wherein the matrix of 
correlation results may be spanned by an index indicating an 
offset from a nominal carrier frequency and also by an index 
indicating code phase, and wherein the selected phasor may be 
the phasor having the maximum magnitude of all the elements of 
the set. Further, the means for providing the matrix of 
correlation results may include means, responsive to a series of 
signal fragments, for performing a coherent integration of each 
of the series of signal fragments, and also means, responsive to 
the coherent integrations, for providing a non-coherent 
integration in which the phasor results of the coherent 
integrations are combined without regard to phase. Still 
further, the means for performing the non- coherent integration 
may multiply each element of a matrix of correlation results 
provided using a coherent integration of a first signal 
fragment, by the complex conjugate of a corresponding element 



for an immediately preceding signal fragment. Also further, in 
some applications, in providing the matrix of correlation 
results as phasor values and in determining the phase of the 
phasor having the maximum magnitude of all the elements of the 
5 matrix, only at most two phasor values may be held in a memory 

device at any instant of time, and of the two phasor values, 
only the phasor value having the larger magnitude may be saved 
in the memory device before calculating a next phasor value. 

In a third aspect of the invention, a system is provided 
10 including: a transmitter for transmitting a signal having a code 

Q component and a carrier component, and a ranging receiver for 
y receiving the signal and for determining information about the 
m carrier frequency of the signal, the ranging receiver 
4;; characterized in that it includes: means, responsive to 
15-Jj approximately carrier-demodulated received signal fragment, for 
providing a set of correlation results indicating information 
|1J about the correlation of the approximately carrier-demodulated 
l: received signal fragments with a replica of the code component 
p and any remaining carrier component, wherein the set is formed 
20^*^^ using different possible offsets from a nominal carrier 

frequency used to approximately carrier-demodulate the received 
signal fragment, and further wherein each element of the set is 
provided as a phasor having a phase and a magnitude; and means, 
responsive to the matrix of phasors, for selecting the phasor 
25 having a magnitude distinguishing it from all the other elements 

of the set, and determining the phase of the selected phasor, 
and for providing information about the carrier frequency based 
on the phase of the selected phasor. 

Further in accord with the third aspect of the invention, 
30 the system provided by the invention may also include a 

computing resource external to the ranging receiver, and the 
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apparatus may communicate information to the computing facility 
via a wireless communication system and the computing facility 
may provide at least some of the computation needed either to 
provide the set of correlation results or to select the phasor. 

5 The invention provides fine frequency estimates of the 

carrier frequency used to convey ranging signals immediately 
after initial acquisition; such estimates can be used for 
tracking or as the initialization of more sensitive dine 
acquisition algorithm, and so simplifying the search for the 
10 carrier frequency. The invention thus allows removing a fine 

acquisition stage from a ranging receiver, saving in some 
O embodiments for example as much as one or two seconds for 
Q acquisition in typical ranging receivers, or simplifying the 
4^ operation of such a stage, reducing the computational burden in 
15,^^. some embodiments for example by as much as a factor of up to 
fifty. 

H Brief Description of the Drawings 

ff The above and other objects, features and advantages of the 

invention will become apparent from a consideration of the 
20 subsequent detailed description presented in connection with 

accompanying drawings, in which: 

Fig. 1 is a flowchart of the invention; 

Fig. 2 is a block diagram/ flow diagram of a GPS receiver 
in which the invention can be implemented; and 

25 Fig. 3 is a block diagram/ flow diagram of an apparatus 

that implements the invention. 
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Best Mode For Carrying Out The Invention 



The invention is described below as a subprocess of the 
acquisition stage of a ranging receiver such as a GPS receiver 
and when specific numerical examples are given, civil GPS 
5 signals (C/A code signals) are used, but it should be understood 

that the invention can be used in any receiver of a spread 
spectrum signal, not necessarily receivers used to provide 
navigation solutions, and the invention is of course not 
restricted to use with civil GPS (C/A code) signals. The 
10^ invention provides a relatively precise estimate of the offset 
n of the carrier frequency from a nominal carrier frequency, the 
H offset attributable to one or another of various causes 
gi including Doppler shifting or clock inaccuracies (the clock of 
'p: either the ranging receiver or the clock of the source of the 
15^11 spread spectrum signal, i.e. a satellite in case of a GPS 

y receiver) . The invention achieves its object by performing, as 
111 part of the processing of the acquisition stage, a modified 

correlation (modified as compared to the correlation performed 
O in prior art ranging receivers) of a received signal fragment 
20" with a replica of the PRN code known to be used by the source of 
the signal . 

As explained above, the acquisition process determines the 
code phase (the relative phase of the replica and the received 
signal fragment) and also the frequency of any complex 
25 sinusoidal modulation remaining after carrier wipeoff assuming a 

nominal carrier frequency, i.e. it also determines a remaining 
frequency offset from the nominal carrier frequency. 

The code phase is determined according to the peak of a 
correlation function. As explained above, the length of the 
30 received signal fragment used in performing the correlation is, 

according to the prior art, restricted typically to 2 0 ms 
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because modulation of the received signal fragment by the 50 Hz 
navigation data corrupts an ordinary correlation using a longer 
signal fragment . 

A process according to the present invention for arriving 
at a more precise estimate of the remaining frequency offset /Sf 
(from a nominal carrier frequency f^) begins with a signal 
fragment after both code wipeoff and nominal carrier frequency 
demodulation by each coarse frequency. The description that 
follows will use equations (1) ~ (6) that are valid only after 
demodulation by the correct code phase and also by the correct 
coarse frequency, i.e. the code phase and coarse frequency 
offset for which a process of coherent and non- coherent 
integrations in performing a correlation of the signal fragment 
with a code replica and a carrier replica (i.e. a sinusoid at a 
frequency equal to the coarse frequency offset) , as described 
below, gives the largest correlation result. For other code 
phase and coarse frequency pairs, more general equations apply, 
equations that reduce to equations (1) - (6) for the correct code 
phase and coarse frequency pair. 

For purposes of the present invention, the signal fragment 
can be modeled as, 

xin) = a„e'-^- +x„,,,(n) = a„e^-W..Ar.)»A, ^^^^^^^ 

where the remaining frequency offset A/ is presented as a 
combination of a coarse frequency offset Af^ with some 
granularity A/^^^^^^ (e.g. 100 Hz) and a fine frequency offset A/^ , 
which is preferably in the range (-A/^^^,^^/2, +A/;^^,^^/2) . The code 

wipe-off in the context of the present invention is to be 

understood as the multiplication of the signal fragment by the 
replica code without integration. If the code phase is correct 
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then the code is wiped off and only modulations due to 
navigation data bits and residual Doppler remain. If the coarse 
frequency offset guess during the acquisition process is 
correct, then the correlation peak will be found and the fine 
5 frequency offset A/} can be estimated. 

In the following description of the invention, the noise 
component x^^^^{n) is omitted for clarity. Moreover, the 
amplitude in eq. (1) should be understood to be the amplitude 
of the signal fragment x(n) , i.e. the amplitude of the signal 
lOy, after code wipeoff using the correct code phase, and after 
y frequency demodulation using a suitably good estimate of the 
H coarse frequency offset Af^ (a good estimate in that it is 
jz substantially the best guess possible using a step size of 50 
^ Hz) . If the estimate of the coarse frequency offset yields a 
15:. frequency that is far from the true carrier frequency (including 
p the Doppler, etc.) or if the code phase is not correct, then the 
H signal value x(n) is small, i.e. it is suppressed. The 
H description that follows does not address such cases. The 

problem addressed by the invention is to find a fine frequency 
20 offset after the acquisition search determines the correct code 

phase (i.e. code wipeoff is assumed) and some coarse frequency 
estimate (so that the signal arriving at a logic module 
implementing the invention is assumed demodulated by the 
estimated coarse frequency) . Of course in actual practice, the 
25 inputs to a logic module that implements the invention will be a 

number of different correlations, each based on a coherent 
integration of the received signal fragment with a particular 
code phase and a particular coarse frequency offset* 

The integration (summation) of an ordinary correlation 
30 calculation subsamples and thereby enhances the signal 
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component, each subsample corresponding to a different increment 
of t4f , i.e. A/} , 24/}, 3 A/} , . . . , ^Mf / where is the number 
of increments of t^f^ included in the ordinary correlation. 

After demodulation by the factor e'^'^^^^ to account for the 
5 coarse frequency offset, each siibsample indicated by the index 

n^n^N ^^n^ , can be expressed as, 

B2=0 

where is the number of samples in the coherent integration 
J:[ (indicated by eq. (2)) and is usually some integral multiple of 
10^^ the PRN code epoch. The result, X^i) / is the (integrated) 
4» signal at a discrete time indicated by ; y{n^ is, more 
r|J specifically, a fragment of the received, code-wiped signal, iV^ 

H samples long, correlated with the sinusoid e"^^^"'^^"""'^^' . (The 
y: index n is used to indicate an instant of time before the 
IBllJ coherent integration, with n incrementing with each new sample, 
p^j; For example, if there are two samples per chip then the sampling 
rate is 2*1,023,000 samples per second, and so n increments 
2*1,023,000 for each second.) After the coherent integration, 
which is an integration over an integral multiple of the code 
20 epoch length (i.e. an interval of length k*1023 chips for some 

integer k) , a second time index is used, one that increments 
every coherent integration length. For example, if the coherent 
integration length is just one code epoch (i.e. 1023 chips), 
then the integrated sampling rate will be 1000 samples per 
25 second; i.e. with the coherent integration, every 1023 samples 

are accumulated into a single, integrated sample. 
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Bear in mind that, as mentioned above, in actual practice 
the quantity y(n^) is calculated for a particular code phase and 
a particular coarse frequency offset. If the code phase is 
indicated by an index p and the coarse frequency is indicated 

by an index m , then y{n^) is more properly written as yp^^^in^) , 
but for simplicity of notation, the p,m dependence is 
suppressed in what follows. 

Multiplying each sample after the coherent processing, 
given by eq. (2) , by the complex conjugate of the nearest 
neighbor sample occurring earlier in time, provides, 

which compensates for the time (code epoch) dependence . 
Defining ^„^^g,^^^^^+^ , and using the well-known Euler's theorem, 
the quantity z(«j) can be written as, 

z(w,) = ^„^/^^^^^^' =J„Jcos(2;r^^^ . (4) 

The amplitude a^^ of the sinusoidal quantity z(kj) in general 

can vazy with time (since it depends on , which indicates an 
instant of time) , due to the BPSK data modulation, while the 
other terms of the quantity z(wj) (the terms in square brackets 
in eq. (4) yielding the phase of the phasor quantity z{n^)) do 
not depend on time. And here is the key to understanding the 
invention: one can find the phase of the phasor quantity z{n^) 
easily enough because that phase does not change with time, and 
since that phase has information about the sought-after fine 
frequency offset A/} , we can use the phase to find the fine 

frequency offset A/^. . 
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Now, according to the invention, a more precise estimate of 
the frequency offset than results from an ordinary correlation 
is provided by what is here called a (supplemental) non-coherent 
processing stage that integrates (sums) the absolute values of 

real part of the signal z(n^) for different times each indicated 

by a different value of i\ , and, separately, the imaginary 
parts, as follows: 

and 

c('-> = ^£|lm(<«„.)| - E I sin(2;rA/?V,A0 1 =| sin(2;rA/ZV,A0 i ""^'l | , (6 ) 

where N^^ is the number of integrated samples in the noncoherent 
integration, i.e. the noncoherent integration length. The value 
used for N^^ depends on the signal strength; it is smaller for a 
higher signal strength, and could even be as small as one. 

As explained above, the amplitude and so the amplitude 

a^^ and likewise the amplitude a^^ are amplitudes for a 
particular code phase p and a particular frequency offset m , 
both p and m being integers and serving as indices, but the 
indices p and m are suppressed in eqs. (1) - (6) above for 

clarity of notation. Thus, the amplitude a^^ is really 

shorthand for a^^^^^^ where the index refers to a particular 
code epoch (and so to a particular instant of time) . Thus also 
z(wj) is shorthand for ^^^^(n^) , and so c^''^ is shorthand for c^^^ and 

similarly for c^^^ . 
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Now according further to the invention, after performing 
the two non-coherent integrations given by eqs. (5) and (6), and 

constructing the matrix c^^^ according to. 



c,^=ci'l+jc^!^^ (7) 

p,m p,m J p,m ^ * 

(now again making express the p^m dependence of c^^m 
eqs. (5) and (6)), the set of all 



is searched to find which is the maximum, indicated here as 



max 

p,m 



PyM 



attains its 



"^^^ indices p and m for which 

maximum indicate the code phase and coarse frequency offset. 
According to the invention, the real component c^^^^ with pytn for 

which Cp J^ attains its maximum, and the corresponding imaginary 

component c^^^ are then used to determine the fine frequency 

offset Isfj- on the basis that these two quantities indicate the 

phase of a phasor c^j^^c^p^^+jc^pj^ (for the best p^m , i.e. for the 

p^m for which |c^,„| attains its maximum) , a phase that depends on 

the fine frequency offset Af^ . 

Thus, at each pair of values of frequency and code phase 
that would be searched in a conventional acquisition, the 
invention provides not one correlation output (including 

coherent and noncoherent part), but two outputs, c^*]^ and c^p^^n - 
The strength of the response is calculated using the sum of the 
squared values of the two outputs, and the peak in the strength 
so calculated corresponds to the correct frequency and code 
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phase pair. But in addition, the two outputs c^p^^ and c^^^^ give 

information useful in determining a fine frequency offset, i.e. 
a frequency offset fine enough to estimate the carrier frequency 
to within better than 50 Hz. (After the coherent integration, we 
5 do have complex correlation values and in principle it is 

possible to obtain a fine frequency offset directly from the 
output of the coherent integration, but in practice, in weak 
signal conditions, it is necessary to essentially amplify the 
correlation result using the non-coherent integration process of 
10 the invention. As mentioned above, a longer- interval coherent 

U integration will not work well because of the modulation of the 

C3 received signal by the navigation data.) 

01 In summary, in a conventional noncoherent integration 

process, inputs x{n) (correlations of signal fragments after 
15j^" code wipeoff and coarse acquisition) are complex numbers a„-^j\ 
(more properly <^n,p,m j\,p,m ' indices p and m, indicating 

^ code phase and frequency offset respectively, are dropped here 
13 for simplicity of notation) , and noncoherent integration is 

performed (in a calculation to amplify the correlation of the 
20 signal fragment with a signal of each of several different code 

phases and each of several different frequency offsets) as an 

integration of absolute values + j^^j or squared values 

\a„-\-jb„f . In the present invention, in contrast, noncoherent 

integration is performed by: first, multiplying each input by 
25 the complex conjugate of its neighbor, i.e. forming the product 

i^n j\) ' i^n-i ' second, noncoherently integrating, 

separately, the imaginary part (^A-j-^A-i) t:he real part 

(^n^«-i +^A-i) • '^^^ result is two values {c^''\c^'^) (actually 
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(4'^m^^p,m) / t)ut the dependence on the code phase index p and 
frequency offset index m are dropped for simplicity of 
notation) for the correlation at each frequency and code phase 
pair of the search. The largest of the sum of the squares of 
c^''^ and c^'^ (or the square root thereof) is used to determine the 
fine frequency offset, based on the phase of the phasor c^''^+jc^'K 
The fine frequency offset so determined, when combined with the 
coarse frequency offset (determined as in the prior art to a 
precision of for example +/-50 Hz) and with the nominal 
frequency, yields a precise value (for example, to within 1-2 
Hz) for the carrier frequency (including Doppler and other 
deviations from the nominal carrier frequency) . 

Eqs. (5) and (6) allow only an ambiguous determination of 
the fine frequency offset A/^ ; for a properly designed 

acquisition system (i.e. having a proper choice of frequency 
spacing), four solutions are possible: 



(9) 



(10) 



and 



(11) 



(12) 



where c*" is shorthand for c^^^ for the p,m for which 




similarly for c 
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The estimate of the fine frequency offset is ambiguous 
(having the four possible results given by eqs. (9) -(12)) 
because of taking absolute values of imaginary and real 

components. Nonetheless, the uncertainty in the fine frequency 
5 offset is greatly reduced using the invention. The ambiguity 

can be resolved by one or another method executed by the 
acquisition stage prior to handoff to the tracking algorithm, or 
the tracking algorithm can be initiated with all four possible 
solutions as input. 

10 Referring now to Fig. 1, a flowchart of the invention is 

shown as beginning with a first step 11 in which a received 

O signal fragment x{n) , at a time indicated by the variable n , is 

m correlated with a replica code at a number of different code 

Tl phases indicated by index p and also correlated with (in 

IS'm effect) a number of different replica carrier sinusoids, each 

y including a different coarse frequency offset indicated by index 

ii 

m , producing an output Vp^jr^in^) for each different code phase and 
[if coarse frequency offset, the index incrementing over every 
H A^^ samples, whereas the original index n increments with each 

20 sample. Next, in a step 12, each y^^S^) multiplied by the 

complex conjugate of its neighbor, i.e. by >'p,;„(«i -1) / producing 
Zp^{n^) having real and imaginary components, for all different 
p.m. Next, in a step 13, the real and imaginary components of 
the quantities for all p^m are each (separately) 

25 noncoherent ly integrated over N^^ values of Wj (each value of «, 

corresponding to A^^ samples), as indicated by eqs. (5) and (6). 
The result is a quantity c^^ having real and imaginary 
components c'p^ and c^'^^ respectively. Next, in a step 14, the 



-18- 



^p,m searched over all p,m to find the p,m for which \Cp.m\ 
{or, equivalently, \cp^n^) is a maximum, yielding a "best" p,m 
and corresponding c^'^^ and c^p^^n- Finally, in a step 15, 
(ambiguous) information about the fine frequency offset A/^ is 

determined based on the ratio c^p^jn^^^P,mf eqs. (9)- 

(12) . As explained above, the ambiguity in the fine frequency 
offset can be resolved by further processing within the 
acquisition module, or the ambiguous information can be provided 
to the tracking module and the ambiguity can be resolved there. 

Still referring to Fig. 1, in essence the method includes a 
first step 100, including substeps 11, 12, and 13, which 
provides a matrix 104 of correlation results indicating 
information about the correlations of the approximately carrier- 
demodulated successively received signal fragments 102 with a 
replica of the code component and a replica of any remaining 
carrier component (i.e. a sinusoid having a frequency equal to a 
coarse offset from the nominal carrier frequency) , wherein the 
matrix is spanned by the index p indicating code phase and the 
index m indicating an offset from a nominal carrier frequency 
used to approximately carrier-demodulate the received signal 
fragments, and further wherein each element of the matrix 104 is 
provided as a phasor (Cp^m) / and a second step 104 of responding 
to the matrix 104 of phasors and determining the phase of the 
phasor having the maximum magnitude of all the elements of the 
matrix, and so providing a more precise value for the carrier 
frequency of the carrier than what is provided by the offset for 
the phasor having the maximum magnitude - 

Referring now to Fig. 2, an exemplary embodiment of a 
ranging receiver 2 0 in which the invention can be implemented is 
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shown receiving a signal s (t) corresponding to a transmitted 
signal St(t), given by, 

St(t) = A C(t) P{t) cos(wt). 

The receiver 20 includes an acquirer 21 (in which logic code 
5 implementing the steps of the invention as indicated in Fig. 1 

is hosted) for acquiring the signal, and a tracking module 24 
for making adjustments necessary to track the acquired signal, 
i.e. the maintain acquisition. The received signal s(t) 
includes a code component C(p(t) at a phase 9 relative to the 
10 phase of the code component of the transmitted signal, a carrier 

Jl! component Axcos (w' t+cpcar) / the cpcar indicating that the carrier 
O component of the received signal is out of phase by an amount 
m (pear with respect to the transmitted signal, and a data component 
P(t), again in general differing in phase from the data 
15.0 component of the transmitted signal although not expressly 
:^ indicated. The carrier frequency w' is different than the 
|\e transmitted frequency w (w' is shifted in frequency from the 
fj; nominal carrier frequency w) because of the relative motion of 
0 the transmitter and receiver, i.e. because of Doppler shifting, 
20^'''' and also because of inaccuracies in the receiver and transmitter 
clocks. A replica of the code component C<p(t) is generated in 
the receiver (in either the acquirer 21 or the tracking module 
24, depending on the stage of operation) with some phase 9. In 
routine receiver operation, after the receiver has acquired a 
25 signal, in order to extract the data component the receiver 

mixes the received signal with a sinusoid at a received carrier 
frequency w' , and also mixes the signal with a replica code 
C(p(t) at a replica code phase cp that puts the replica code in 
phase with the transmitted code. (The phases are different 
30 because of the transmission time from the satellite to the 
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receiver and also because the clocks in the receiver and 
transmitter are not synchronized.) 

The received carrier frequency w' and the correct replica 
code phase 9 are determined by an acquirer 21. To acquire the 
5 signal s(t), and so to determine the received carrier frequency 

w' and replica code phase 9, the receiver provides a portion of 
the received signal s{t) to a radio- frequency (RF) to baseband 
module that includes an RF downconverter and an analog-to- 
digital converter (ADC) , providing an approximately baseband 
10 signal (n) (the index n indicating a sample made at a 

particular time and so indicating the particular time) , which 
y still includes the code component C (n) , The RF to baseband 
M module provides in-phase and quadrature versions of the 
^ received, downconverted signal, and provides the two versions to 
15|1 respective ADC modules. One ADC module then provides an in- 
"'^ phase digitized, downconverted received signal s'i(n), and the 
H other provides a quadrature digitized, downconverted signal 
s'gCn) . The combination of two such outputs is treated as a 
ry complex entity in what follows. 

Instead of an RF to baseband down- converting followed by an 
analog to digital conversion, it is also possible to use an RF 
to intermediate frequency conversion followed by a conversion to 
baseband combined with ADC. 

Once the acquirer 21 precisely determines the carrier 
25 frequency w' , including any shifting 6w away from the nominal 

carrier frequency w, and also the replica code phase (p, the 
receiver 20 can extract the data signal P(t) from the received 
signal s(t) by for example first mixing the received signal s(t) 
with a sinusoid at the precisely determined carrier frequency 
30 w' , then mixing the resulting signal with a synchronized (in 
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phase) replica code signal C(p(t), and finally performing an 
integrate and dump, leaving only the data signal P(t). 
Sometimes the data extraction is performed within the tracking 
module 24, as is indicated in the embodiment shown in Fig, 2. 
5 The tracking module then not only extracts the data signal P(t) 

but also examines it to determine whether to make slight, 
tracking adjustments to the carrier frequency w' or the replica 
code phase cp so as to keep the receiver 2 0 tuned to the signal 
s(t) . 

10 Still referring to Fig. 2, the signal s' (n) , referring to 

JJ; both s'i(n) and s'Q(n), includes both a code component C(n) as 
l3 well as a residual carrier component depending on the shift 6w. 
g'l As explained above, the shift 6w in some applications is 
tr. typically anywhere in a bandwidth of 12 kHz about the nominal 
15|| carrier frequency w. The code C(n) is typically a pseudorandom 
iT^ sequence of bits of a predetermined length. In case of civil 
fli GPS applications, a code-length of 1023 bits is used. Since the 

code does not itself convey information and is instead merely a 
D device by which the bandwidth of an information-bearing signal 
20' is spread over a much wider bandwidth in a way that allows other 

information-bearing signals to also be spread over and use the 
same bandwidth, a bit of such a code is referred to as a chip, 
instead of a bit. The chip rate is typically much greater than 
the bit rate of the information-bearing signal P (n) . In GPS 
25 applications, the chip rate is typically 1023 chips per 

millisecond, whereas the data P(n) is nominally 50 bits per 
second . 

Still referring to Fig. 2, the acquirer 21 is shown as 
including a coarse carrier wipeoff module (CCWM) 22 including a 
3 0 numerically controlled oscillator (NCO) , a processing module 23 

(hosting the code logic for practicing the invention) , a control 
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module, a scanning module, and switching modules used to switch 
from acquiring to tracking after a signal is satisfactorily 
acquired. The control module operates the switches SWl and SW2 
so that once the signal has been acquired (with respect to both 
5 carrier frequency and code phase) , the acquisition module is 

switched out of the circuit including the NCO, and the tracking 
module 24 is switched into that circuit. 

The approximate baseband signal s' (n) is provided to the 
CCWM 22, which mixes s' (n) with a complex sinusoid (i.e. an in- 
10 phase and quadrature sinusoid are mixed with the in-phase and 

quadrature approximate baseband signals s'i(n) and s'Q(n)), at 
0 one of a number of several distinct frequencies in a set of 
p coarse frequency offsets spanning a range of frequencies about 
0'; the nominal frequency encompassing all possible frequency 
isjj shifts. For example, in the case of a GPS receiver, the CCWM 22 
^ would first mix the input s' (n) with a signal at 6 kHz below the 
i:^ nominal carrier frequency, to provide x(n) . Later, at the 

command of the analysis module 23d, the CCWM 22 would mix the 
fij same signal s' (n) at a next coarse frequency, 5 kHz below the 
20H nominal carrier frequency, and so on, each time producing a 
different signal (complex) x(n) including a different coarse 
frequency component. 

As will be explained below and as mentioned above, the 
processing module 23 includes components for performing the 
25 steps indicated in Fig. 1. The output of the processing module, 

which is used as the input to the tracking module 24, is 
preferably an unambiguously determined carrier frequency w' and 

a code phase 9. 

Note that according to the invention, a carrier wipeoff 
30 using each of the coarse frequencies is needed by the processing 

module 23 before the processing module can begin the analysis 
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used to acquire a signal. Thus, the processing module does not 
pass control to the tracking module 24 until after having the 
CCWM 22 perform an approximate carrier wipeoff with each coarse 
frequency spanning a range of frequencies presumed to include 
the actual to-be-determined carrier frequency, and then 
performing the analysis used to indicate a precise value for the 
carrier frequency according to the steps indicated in Fig. 1. 
In the course of performing the analysis, as explained above, 
the phasors c^^^ corresponding to the different coarse 
frequencies are constructed and that having the largest 
magnitude is determined. To find the c^^^ for which \cpj^ is a 

maximum, either each c^^ can be saved and the resulting array 
of searched after all c„^ are constructed, or instead only 

p,m p,m ' 

the for which \c^A is so far the maximum can be saved. 

Referring now to Fig. 3 and also to Fig. 1, the processing 
module 2 3 of Fig. 2 is shown in more detail as including modules 
implementing the steps indicated in Fig. 1. A module 31 (Fig. 
3) performs the coherent integration of step 11 (Fig. 1), 
storing the results of each coherent integration for a 
particular code phase (indicated by an index p) and coarse 
frequency offset (indicated by an index m) in a memory device 
31a. When all of the coherent calculations are completed (i.e. 
for all p and all m, the range in p and m sufficient to span all 
possible code phases and all possible carrier frequency shifts) , 
a neighbor multiplication module 32 then performs the neighbor 
multiplication step 12. Next, a non-coherent integration module 
33 performs the non-coherent integration step 13, amplifying the 
results of the coherent integration. A single matrix of results 
(each results being a complex number) is provided as an output 
of the non-coherent integration module 33; the matrix is spanned 
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by the code phase index p and the coarse carrier frequency index 
m. A maximum finder module 34 then searches the matrix and 
finds the p and m for which the magnitude of the complex matrix 
elements is a maximum (or, equivalent ly, the square of the 
5 magnitude) . Finally, a solver module 34 determines the phase of 

the phasor corresponding to the complex elements for the found p 
and m, extracts information about the fine frequency offset 
(according to step 14) , and, in the preferred embodiment, 
performs further processing to eliminate the ambiguity of the 
10 results described in connection with eqs. (9) - (12) . The output 

l^s. of the solver module is therefore, in the preferred embodiment, 
y the precisely determined carrier frequency (including both the 
M coarse and fine frequency offsets) and also the code phase, both 
of which are then used by the tracking phase of the receiver. 

15^:^ Still referring to Fig. 3, from another perspective, the 

■« invention is in essence an apparatus including: a first module 
.H= 3 00, responsive to successive approximately carrier-demodulated 
H received signal fragments 302, for providing a matrix 304 of 

correlation results indicating information about the correlation 

2 OH' of the approximately carrier-demodulated received signal 

fragment with a replica of the code component and a replica of 
any remaining carrier component, wherein the matrix is spanned 
by an index p indicating code phase and an index m indicating a 
first offset from a nominal carrier frequency used to 

25 approximately carrier-demodulate the received signal fragment, 
and further wherein each element of the matrix is provided as a 
phasor Cp,m having a phase and a magnitude; and a second module 
3 06, responsive to the matrix of phasors 3 04, for determining 
the phase of the phasor having the maximum magnitude of all the 

30 elements of the matrix, and for providing information about the 

carrier frequency based on the phase of the phasor having the 
maximum magnitude of all the elements of the matrix. In the 
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preferred embodiment, the first module 300 includes: a first 
sub-module 31, responsive to a series of signal fragments, for 
performing a coherent integration of each of the series of 
signal fragments, preferably by multiplying each element of a 
5 matrix of correlation results provided using a coherent 

integration of a first signal fragment, by the complex conjugate 
of a corresponding element for an immediately preceding signal 
fragment; and also a second sub-module 32, responsive to the 
coherent integrations, for providing a non-coherent integration 
10 in which the phasor results of the coherent integrations are 

^ combined without regard to phase. 

It should be understood that the in performing the neighbor 
multiplications for forming the noncoherent sums, although the 
% invention has been described above in which neighbor 
IsJi multiplications are multiplications of the first and second 

J'' coherent integration result, then the second and third, then the 

third and fourth and so on, i.e. (1x2), (2x3) , (3x4) , the 
H invention also comprehends using other neighbor multiplications, 
g such as (1x2), (3x4), (5x6), ... , or (1x3), (2x4), (3x5), (4x6),... 
2rf"'' (although the first type of neighbor multiplication (i.e. (1x2), 
(2x3) , (3x4) , ... ) is preferred by the inventor. 

It should also be understood that although the above 
description has been in terms of constructing a matrix c^^^ of 
correlations for p,m pairs (the matrix given by eq. (7) , with 

25 values for the real and imaginary components of each element 

given by eqs. (5) and (6) respectively after demodulation of a 
signal fragment using the correct code phase and coarse 
frequency as explained above) , the invention is not limited to 
implementations in which all of the matrix elements are held in 

30 a memory device at any one time. In particular, it is usually 
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preferable to perform all the integrations for a first p,m 
pair, then for each next such pair, perform the calculations and 
compare the results to those for the previous pair and keep only 
the largest result (of course also keeping note of the p^m pair 
5 that gave rise to the current largest result) . 

Scope of the Invention 

It is to be understood that the above-described 
arrangements are only illustrative of the application of the 
principles of the present invention. In particular, it should 
10 be noted that the invention is in no way limited to use in a 

h ranging receiver as illustrated in Fig. 2, which is provided 
O only as an example of a ranging receiver in which the invention 
m can be implemented; as another example, besides using a NCO to 
=f;= Wipe off (approximately demodulate) the carrier (by removing a 
15j| coarse frequency) , discrete Fourier transforms (DFTs) are used. 
:\ Numerous modifications and alternative arrangements may be 
III devised by those skilled in the art without departing from the 
iT. spirit and scope of the present invention, and the appended 

ty 

p claims are intended to cover such modifications and 
2 (f^ arrangement s . 
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